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A benchmark problem for numerical relativity has been the head-on collision of two black holes 



starting from the "Misner initial data," a closed form momentarily stationary solution to the con- 
straint equations with an adjustable closeness parameter fiQ. We show here how an eclectic mixture 
of approximation methods can provide both an efficient means of determining the time development 
of the initial data and a good understanding of the physics of the problem. When the Misner data is 
chosen to correspond to holes initially very close together, a common horizon surrounds both holes 
and the geometry exterior to the horizon can be treated as a non-spherical perturbation of a single 
Schwarzschild hole. When the holes are initially well separated the problem can be treated with a 
different approximation scheme, "the particle-membrane method." For all initial separations, nu- 
merical relativity is in principle applicable, but is costly and of uncertain accuracy. We present here 
a comparison of the different approaches. We compare waveforms, for 1 = 1 and I = 4 radiation, for 
different values of /io, from the three different approaches to the problem. 
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I. INTRODUCTION 

The advent of LIGO and VIRGO §, and other grav- 
itational radiation detectors in the coming years, adds 
motivation to achieve a better understanding of the 
physics of black hole collisions, physics that is already 
interesting for many other reasons. Not the least of these 
is that black hole collisions were the earliest testbed for 
numerical relativity, the solution of the differential equa- 
tions of Einstein's theory with numerical codes. Smarr 
and Eppley f|J|f|, more than 15 years ago, computed 
the radiation waveforms for the axisymmetric problem of 
two holes, starting from rest and falling into each other 
in a head-on collision. Partly to gauge the progress of the 
past 15 years, both in computing machinery and in its ap- 
plication to the problem, the problem has recently been 
reconsidered and recomputed H||. One of the reasons 
for the importance of this problem is that it is a starting 
point for a major initiative in state-of-the-art computing, 
the Grand Challenge project @ to compute the spiraling 
coalescence of binary black hole systems in the absence of 
symmetries. This project will rely on both large scale nu- 
merical computations to solve the fully nonlinear Einstein 
equations and on a new set of semi-analytic techniques, 
such as those presented here, to gain a physical under- 
standing of the numerical results and to provide testbed 
calculations to verify the new generation of numerical 
codes under development ||. 

The technical difficulties of solving the dynamics of 
black hole interactions are daunting, as will be clear in 
the limits of progress reported below. The difficulties 
make it important, or even necessary, that the numerical 
work be supported by other approaches. In particular, 
the numerical challenge would be greatly aided by ana- 
lytic or semi-analytic methods of computing black hole 
solutions in cases where such methods can be devised, so 
that those developing numerical codes would be guided. 
Most important, distinctly different methods of finding 
solutions can provide a measure of the errors, especially 
the systematic errors, in the numerical schemes. 

There is a very separate motivation for developing al- 
ternative methods of finding solutions. While numerical 
treatment gives us the data representing the physical pro- 
cess, analytical methods provide a structure for under- 
standing these data, for perceiving what is interesting, 
what is expected and is new, and what generalizations 
are plausible, and what questions should be asked next. 

We report here on work done with just such motiva- 
tions. The problem is the head-on collision of two holes, 
the classic problem of Smarr and Eppley. The initial 
data is the analytical solution of the Einstein initial value 
equations given by Misner This solution represents 
two symmetric, momentarily stationary, "throats" with 
a proper separation L, in a spacetime with ADM mass M 
(thus representing two momentarily stationary holes each 
of mass roughly M/2) [jioj. In this solution the dimen- 
sionless measure of initial separation L/M of the throats 



is given by a parameter /io. 

Initial data with large values of fio represent infall sce- 
narios starting from large separations. Small [Iq cor- 
responds to very different collapses starting with close 
throats. Three distinctly different methods are used to 
study the problem: (i) Over a broad range of separations, 
the more-or-less established approach to numerical rela- 
tivity, the numerical solution of the Einstein equations in 
3+1 form, is used to find the spacetime evolved from the 
initial data, (ii) For large separations ("far limit") we 
use the "particle-membrane" method, which starts with 
a point particle infall treated using perturbation theory, 
and then introduces factors describing the internal dy- 
namics of the black holes, using the black hole membrane 
paradigm, (iii) In the "close limit," when initial separa- 
tions are small, we exploit the fact that the initial ge- 
ometry is nearly spherical outside the initial horizon; the 
evolved spacetime outside the horizon is therefore nearly 
spherical and its development can be approximated with 
the theory of non-spherical perturbations of a single black 
hole. 

The range of validity of the three methods is shown be- 
low and their strengths and weaknesses discussed. The 
conclusion that emerges is that the eclectic approach to 
this problem, using a mixture of distinctly different com- 
putational methods, gives a very robust set of answers, 
with good limits on errors, and much improved under- 
standing of the meaning of the answers. We argue that 
similar approaches should be developed, wherever pos- 
sible, as the forefront of the Grand Challenge initiative 
moves forward. 

The remainder of the paper is organized as follows: In 
Sec. II an outline is given of how the problem is translated 
into a problem suitable for numerical computation, and 
how information is extracted about the outgoing radia- 
tion. The "particle-membrane" method suitable for large 
initial separations is discussed in Sec. Ill and numerical 
results from this method are presented. In Sec. IVA an 
introduction is given to the "close limit," the techniques 
of perturbation theory applied to collisions with small 
initial separations. In Sec. IVB numerical results of the 
close limit are shown and are compared with the meth- 
ods of Sec. II and III. Conclusions are briefly presented 
in Sec. V. 



II. NUMERICAL RELATIVITY FOR THE 
HEAD-ON BLACK HOLE COLLISION 

We use the standard 3 + 1 (ADM) formalism as the 
framework for building a numerical code to solve the fully 
nonlinear Einstein equations, and evolve the Misner ini- 
tial data sets without making any approximations. The 
Cadez coordinate system p"T[ is used, for the most part, 
as it provides natural spherical boundaries at the black 
hole throats and in the asymptotic wave zone where radi- 
ation is extracted. It also utilizes a logarithmic "radial" 
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coordinate to extend grid coverage beyond where radia- 
tion can propagate within a typical run time. However, 
this coordinate system has a singular saddle point mid- 
way between the two holes, making numerical evolution 
quite difficult. For this reason an overlapping cylindri- 
cal coordinate system is used as a coordinate patch to 
cover this region, independently of the Cadez grid except 
at the interface boundaries. The nonsingular cylindri- 
cal metric and extrinsic curvature components are then 
used to correct the corresponding singular Cadez compo- 
nents in the overlapping patched regions throughout the 
evolution. We mention this detail because although this 
procedure was very effective in suppressing numerical in- 
stabilities that can develop at the singular point, it also 
has the effect of introducing low amplitude signals in the 
evolution that have a bearing on the interpretation of the 
gravitational radiation waveforms presented below. This 
work has been discussed extensively in Refs. [|6], |l2| , |l3t , 
where complete details of our numerical calculations and 
results can be found. Here wc focus the discussion on 
the extraction of waveforms from the numerically gener- 
ated spacetime metric. In later sections of this paper, 
we compare the waveforms extracted to those obtained 
using the semi-analytic approaches. 

Our waveform calculations are based on the gauge in- 
variant extraction technique dev elop ed by Abrahams and 
Evans O] and applied in Ref. |l5f to black hole space- 
times. This work, in turn, is derived from the gauge- 
invariant formalism developed by Moncrief . The ba- 
sic idea is to split the numerical spacetime metric into 
a spherically symmetric background and a small non- 
spherical perturbation. First, we expand the metric per- 
turbation in m = spherical harmonics Yiq(8) and their 
tensor generalizations. Then, as described in Appendix A 
of Ref. p5| , Moncrief 's perturbation functions H 2 , G, and 
so forth, are projected out of the numerically computed 
metric components using the orthogonality of the Ye (0) 's 
by performing angular integrals. For example, for each 
^-mode, if the background is written in Schwarzschild co- 
ordinates, the perturbation function H 2 can be computed 
by performing an integral of the numerically evolved ra- 
dial metric function g rr over a 2-sphere 
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The other perturbation functions are computed in a sim- 
ilar way. These perturbation functions will be discussed 
in more detail below in section IV where they are com- 
puted in a different way using analytic approximation 
techniques. 

Once these perturbation functions are known numeri- 
cally, they are used to construct the gauge invariant Zer- 
illi function ip. Following Moncrief |l(| we define the 
following quantities 



fci = K + SrG,. 
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where 5 = 1 — 2M/r. Then the quantity defined by 
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is gauge invariant and satisfies the Zerilli equation de- 
scribing the propagation of gravitational waves on a black 
hole bac kgro und. With our choice of normalization given 
by Eq. (2.4), the asymptotic energy flux carried by the 
gravitational wave is given by 



dE__J_ (d± 
dt ~ 32^ \dt 



(2.6) 



for each £ mode. 

It is important to note that this procedure for extract- 
ing the Zerilli function ip differs slightly from the pro- 
cedure used in Sec. Here we imagine the spacetime 
at a large distance to be well represented by spherical 
plus non-spherical pieces, corresponding to various £— 
poles. Although we assume such nonspherical pieces to 
be small enough that nonlinear terms in the non-spherical 
pieces can be neglected with impunity, the procedure de- 
scribed above simply lumps "everything not spherical" 
into the perturbation terms, whether it is small or not. 
In this sense our numerically extracted waveforms could 
be considered "nonlinear." In the semi-analytic approx- 
imation methods discussed below higher-order terms are 
explicitly dropped everywhere in the spacetime, so that 
only the "true" first-order part is included. Of course 
we expect the methods to agree in the regime where the 
approximations we use are valid, such expectations are 
borne out as we show in section IV B. 

We have extracted both the £ = 2 and £ — 4 wave- 
forms, using the method described above, at radii of 15, 
20, 25, 30, and 35 M, where M is the ADM mass of the 
spacetime. By comparing results at each of these radii we 
are able to check the propagation of waves and the con- 
sistency of our energy calculations. Possible sources of 
numerical errors in the gravitational wave signals include: 
(i) truncation error, (ii) artificial diffusion, (iii) grid spac- 
ing that increases exponentially with "radius" , (iv) severe 
pathological behavior arising from the singularity avoid- 
ing (maximal) time slicing, and (v) coordinate patch im- 
plementation. These effects have been discussed in detail 
in Ref. |l3| ] where we have performed a number of con- 
vergence studies and demonstrated the robustness of the 
£ = 2 waveform extraction to a few percent under signif- 
icant changes in computational parameters such as grid 
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resolution, patch width, numerical diffusion, etc. The 
most prominent numerically induced feature in the £ — 2 
waveform is a slight broadening of the wavelength at late 
times. The £ = 4 waveform is less certain than 1 = 2 
(especially at low values of /io) due to its much smaller 
amplitude which makes it difficult to extract from the 
background noise level. Moreover, the more complicated 
angular distribution of the £ = 4 component requires 
more angular zones to be fully resolved than the t = 2 
component. Our convergence studies indicate that al- 
though we are unable to determine with great certainty 
the absolute amplitude of the £ = 4 signal, the damping 
rate, frequency and phase of these waveforms are signif- 
icantly more reliable. We shall see the confirmation of 
these points, together with new insights into the numer- 
ical data, from the comparison with the results from the 
semi- analytic approaches. 

III. THE PARTICLE-MEMBRANE APPROACH 

In this section we consider the case in which the two 
black holes are initially well separated — the "far limit.'' 
In this case the evolution of the system can be separated 
into two phases. The first is the infall phase, during 
which the radiation emitted depends almost uniquely on 
the motion of the two bodies, and almost not at all on 
their nature (that is, whether or not they are black holes). 
The second is the interaction phase, during which the 
radiation emitted can depend significantly on the nature 
of the two bodies; it is during this phase that most of the 
gravitational- wave energy is produced. 

We describe a semi-analytic approach, first presented 
in Ref . || , for calculating the gravitational radiation pro- 
duced in the head-on collisions described by the Misner 
data. The approach consists of two parts. In the first 
(particle) part, we calculate the radiation by assuming 
that only one body is a black hole, and that the other is 
a much less massive (structureless) point particle; then 
we extrapolate the results to allow for a mass ratio of 
unity. The first part of the calculation incorporates com- 
pletely the infall phase of the collision, but only partially 
the interaction phase. In the second (membrane) part of 
the calculation, we promote the point particle to a black 
hole, and compute corrections due to the internal dynam- 
ics of that black hole. These corrections are important 
only during the interaction phase of the evolution. The 
first part of the calculation is carried out using black- hole 
perturbation theory; the second is carried out using the 
membrane paradigm for black-hole dynamics. 

The particle-membrane approach was used in Ref. 0] 
to estimate the energy carried off by gravitational radia- 
tion. The result was found to be in excellent agreement 
with the numerical calculation; see Fig. in Sec. Ill A. 
Here we push the comparison one step further: We use 
the particle-membrane approach to calculate the wave- 
forms, which we then compare to those obtained numer- 



ically. 

In Sec. Ill A (infall from infinity) we construct the 
waveform by considering a particle falling in toward the 
black hole from infinity. We shall see that although at 
late times the agreement is very good, the waveform dif- 
fers from the numerical one by some low frequency com- 
ponents at earlier times. This is to be expected, because 
the physical situations considered are very different from 
each other at early times: infall from infinity versus infall 
from a finite distance. These differences in waveforms are 
nevertheless compatible with a good agreement in the to- 
tal energy outputs, as the early part of the infall does not 
generate much radiation. 

In Sec. Ill B (time-symmetric trajectory) we compute 
improved waveforms based on particle motions bearing a 
closer resemblance to the numerical set-up. We consider 
a particle emerging out of the white-hole (past) horizon, 
coming to a stop at some finite distance away from the 
black hole, and then falling back. The waveforms ob- 
tained are compared to the numerically generated ones 
in detail. 



A. Infall from Infinity 

We begin by outlining the particle-membrane calcula- 
tion presented in Ref. |fj. The emphasis will be put on 
those steps which are relevant for the computation of the 
waveforms. 

The starting point is the twenty-three-year old com- 
putation (lTj of the energy radiated by a particle falling 
into a black hole from infinity. The particle has mass mi, 
while the mass of the (non-rotating) black hole is denoted 
77i2. Under the assumption mi <C 7772, the total energy 
radiated can be calculated using black-hole perturbation 
theory p7[ ]. The result is 

E = Km 1 2 /m 2 , k = 0.0104. (3.1) 

To compare this with the numerical results, we must first 
(i) extrapolate to a mass ratio of unity, (ii) take into ac- 
count the fact that in the numerical simulations, the in- 
fall proceeds not from infinity but from a finite distance, 
and (iii) correct for the internal dynamics of the infalling 
black hole. 

To extrapolate to the case mi = 7712, we begin by 
first understanding the origin of the factor mi 2 /"i2 in 
Eq. ([O]). This, we do by invoking the quadrupole for- 
mula for gravitational radiation. If / denotes a typical 
component of the system's quadrupole moment, then 

I-mir 2 , (3.2) 

where r is the distance to the origin (where 7772 is assumed 
to sit). The gravitational-wave luminosity L = dE/dt is 
then given by 
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To obtain the total energy radiated, we integrate this 
over time. The dominant contribution comes from the 
strong-field region, so 



^ ~ ^strong field ^strong field ' 



(3.4) 



where Ai g j- ron g f^ e ^ is the time required by the particle 
to cross the strong-field region, which has an extension 
r strong field ~ m ^ 80 that A *strong field ~ m 2- Putting 
these into Ecl ( |3~3| ) and Eq. ( |3~4| ) yields E ~ mi 2 /m 2 , 
as in Eq. ( |3.l| ). In particular, we point out that E is 
inversely proportional to 7712, the length scale associated 
with the strong-field region. 

We now repeat this argument having in mind an 
infalling particle which has a mass mi not negligible 
compared with the mass m 2 of the black hole. The 
quadrupole moment must now be replaced by I ~ jir 2 , 
where /i = m\m>xl{rrh\ + m 2 ) is the reduced mass and 
r the relative separation. On the other hand, quite 
independent of the value of mi, the extension of the 
strong- field region as seen by the infalling particle is stil l 

Strong field ~ Ai strong field ~ m ^ Hence Ec l- & 
now leads to E ~ /i 2 /m,2. Finally, taking the limiting 
case mi = mi , and inserting the numerical factor k from 



Eq. (3.1), we arrive at 



E = 2kh 2 /M = kM/8, 



(3.5) 



where M = mi + m 2 is the total mass of the system. 
Notice that, say, if one replaces mi by /z, but m.2 by M 
in (3.1), the result will be different by a factor of 2. We 
note that while there is no definite way to carry out the 
extrapolation, to the extent that a factor of 2 is meaning- 
ful in such arguments, our way of doing it is physically 
motivated. We discuss this point in detail here as in a 
different consideration below, mi will be extrapolated to 
M instead of M/2. 

Next, we modify Eq. ( |3.5| ) to reflect the fact that in the 
numerical simulation, the infall proceeds from a finite 
initial separation (denoted ro) instead of infinity. This 
modification was described in detail in Ref. Q, and we 
will not repeat this discussion here. Suffice it to say that 
an infall from a finite distance produces less radiation 
than an infall from infinity. This is because (i) there is 
less time for the system to radiate, and (ii) the infall 
velocity is smaller at a given separation. The reduction 
in the total energy outp ut c an be obtained by inserting 
a factor J- ro < 1 in Eq. (|3.5|), such that 



E = 2Kf ro iJ 2 /M. 



(3.6) 



An expression for J- ro can be found in Ref. |6J. 

The last step consists of invoking the membrane 
paradigm to calculate the corrections due to the inter- 
nal dynamics of the infalling object, which is now taken 
to be another black hole instead of a structureless point 
particle. These corrections can also be put into the form 
of additional factors to be inserted into the expression for 



E. These factors are all of order unity, implying that the 
internal dynamics of the black hole is unimportant in the 
infalling process, so far as gravitational wave generation 
is concerned, a black hole falling into another black hole 
is not much different from that of a point particle falling 
into a black hole. It is found |J that the largest effect 
comes from the tidal heating of the black hole horizon, 
which is distorted by the gravitational field of the other 
hole as they approach each other. The corresponding fac- 
tor Th is found to be a function of ro, with its smallest 
value ~ 0.86 realized when r — > oo. For the range of r 
covered by the numerical evolutions described in Sec. II, 
the total energy 



E = 2KT ro F h n 2 /M, 



(3.7) 



is reduced just by a few percent by Th- (There are other 
factors with less than one percent effects. They are not 
considered in this paper. For details see Ref. B) 

We compare the predictions of the particle-membrane 
approach with the numerical results in Fig. |l| In the fig- 
ure, the total energy radiated is plotted as a function of 
both L, the proper initial separation (a function of ro), 
and the Misner parameter fiQ. The clustered symbols 
show the energy extracted at various detector locations, 
as obtained from the numerical results. The procedure 
for carrying out this calculatio n wa s described in Sec. II. 
The solid line is a plot of Eq. (3.7). We see that for suf- 
ficiently large init ial s eparations, the agreement is quite 
good. Equation (3.7) overestimates the energy output 
when the initial separation is such that /i < 1-8. This is 
to be expected, because for such small initial separations, 
the two holes are surrounded by a common horizon [ p"8"[ . 
For such cases the particle-membrane approach, which 
is based on an object falling into a black hole, is clearly 
inappropriate. Much better suited is the close-limit ap- 
proach, to be discussed in Sec. IV. 

We now turn to a discussion of the waveforms. In 
Fig. ^| we represent by a solid line the numerically- 
obtained waveform V'numWi as measured by a detector 
at r = 35M, resulting from a collision with Misner pa- 
rameter /j,q — 2.7 (corresponding to L = 6.35M). Wc 
wish to compare this with a waveform obtained using 
the particle-membrane approach. 

To construct a waveform with the particle-membrane 
approach, we begin with ipoo, the waveform resulting 
from the infall of a particle of mass mi into a black hole 
of mass m-2 3> mi, with the infall starting at infinity. The 
waveform is measured by a detector situated near future 
null infinity; it is a function of retarded time t — r* , 



-0OO = i>o 



m 2 



(3.8) 



where r* = r + 2m2 ln(r/2m,2 — 1). The function tfj^ has 
been calculated, using perturbation theory, by several 
authors 

To convert "0oo into a "perturbation-membrane" wave- 
form, ippm, we proceed in three steps. First, we renor- 
malize the amplitude of ipoc in such a way that the total 
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FIG. 1. Comparison of the particle-membrane approach 
computation of the energy vs. the full numerical results. The 
solid line gives the results of the particle-membrane approach, 
while the clusters of symbols represent energies computed nu- 
merically for a series of simulations carried out for different 
values of no- For each value of /xo, the different symbols corre- 
spond to energies extracted at different radii. For small initial 
separations with fj,o < 1.8, with the two holes surrounded by 
a common horizon at the initial time, the particle-membrane 
approach overestimates the energy output. 
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energy carried off by that wave agrees with Eq. (3.7). Sec- 
ond, we replace m 2 in Eq. ( |3.S|) by M, the total mass of 
the spacetime. This is justified by noting that the waves 
reaching an observer situated at a large distance from 
the system are scattered (except for some high frequency 
components) by the curvature of the spacetime whose 
length scale is determined mainly by the total mass M of 
the system. [This is different from what was used earlier 
for the total energy, which depends on the length scale of 
the strong field region experienced by the particle falling 
into the black hole. The generation of radiation energy 
(which strongly depends on near zone physics) involves 
a different scale compared to that of wave propagation 
through the field of the combined holes.] In particular, 
at late times the waveform must be dominated by quasi- 
normal ringing of the final black hole, which has a mass 
mi + TO2 = M. Third, to compare with the numerical 
result f/'num which is calculated for a detector at radius 
r = 35M (instead of at future null infinity where tp p is 
given), we have chosen the phase of ^ pm so as to match 
that of f/'num at late times. The final result is displayed 
as a dotted line in Fig. ^. 

The numerical and particle-membrane waveforms are 
in reasonably good agreement when their respective be- 
havior is dominated by quasinormal ringing. We notice 
that at late times, the wavelength of ■0num becomes longer 
than that of ip pm . This is an artifact of the numerical 
calculation, as was already noted in Ref. (The be- 
havior of i/Win at late times is sensitive to the choice 
of various parameters, such as resolution, used in the nu- 
merical computation.) There is also a significant discrep- 
ancy at early times: Vpm displays a long downward ramp 
corresponding to the slow infall of the particle at early 
times. Because in the Misner data the black holes collide 
from a finite initial separation, this ramp is not found 
in the numerically-obtained waveform. Indeed, the nu- 
merical simulation can cover only cases with relatively 
small initial separation, namely, less than 10 M. For 
larger separations, there are difficulties associated with 
the coordinate system used in the numerical calculation, 
see the previous section. We note that this low frequency 
downward ramp part of i/' 00 does not correspond to much 
radiation energy, as the luminosity is related to the time 
derivative of ip. We therefore have the good agreement 
in the total energy, as shown in Fig. pi despite the dif- 
ferences in the waveforms. 

This discrepancy in the early part of the waveforms 
suggests a refined treatment of the particle part of the 
particle-membrane approach, to which we now turn. 



FIG. 2. Comparison of the particle-membrane approach 
computation of the waveform (dotted line) vs. the full nu- 
merical results (solid line). The waveforms were extracted at 
r = 35M and the initial data are for two black holes separated 
by a distance L — 6.35M (/io = 2.7). 



B. Time-Symmetric Trajectory 

The Misner data describe two Einstein-Rosen throats 
which first fly apart from one another, then turn around 
and fall back in a time symmetric manner. The geom- 
etry is time symmetric about the t — slice at which 
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FIG. 3. A time symmetric particle trajectory used in 
the membrane-particle approximation model. The particle 
emerges from the past horizon at the point E and falls back 
in to the future horizon at a later time. The times labeled r w 
and t s axe the times at which signals reach the detector from 
the events when the particle leaves the white hole and crosses 
the point of time symmetry, respectively. 

the Misner data is given [although the evolution of the 
black hole horizons are not time symmetric (Ref. pl|)]. 
In this subsection we refine our particle-membrane ap- 
proach by considering, instead of an infall from infin- 
ity, a time-symmetric trajectory such as depicted in Fig. 
|3|. In this refined model, the particle emerges from the 
white-hole (past) horizon, travels up to a radius ro, and 
then falls back. The trajectory is a geodesic of the 
Schwarzschild spacctime, and is time-symmetric about 
the turning point r = r<j. 

This problem not only provides a better model for the 
collision of two black holes, it is also interesting in its own 
right. We will see that the particle's emergence from the 
past horizon causes the excitation of the white-hole quasi- 
normal modes. When the particle falls back, it is the 
black-hole quasinormal modes which are now excited, as 
we have seen previously. The white-hole and black-hole 
modes are generated with different amplitudes: Although 
both the background spacetime and the trajectory are 
time-symmetric about t = 0, the radiation pattern is not. 
This should not be surprising, because the radiation is 
calculated using retarded integrals, which break the time 
symmetry. To the best of our knowledge, perturbation- 
theory calculations of time-symmetric trajectories have 
never been carried out before, and the quasinormal ring- 
ing of white holes has never been observed. 

To treat this problem we use Teukolsky's perturba- 
tion formalism [^2), as presented in detail in Ref. f p3| . 
Schematically, the equation to be solved is of the form 



(3.9) 



where 'F represents the perturbation, □ a certain wave 
operator, and p the source term (which is constructed 
from the particle's stress-energy tensor). The solution to 
Eq. (3J) can be expressed as 



* = □: 



(3.10) 



with D" 1 designating the retarded integral, which con- 
sists of three components: (i) the retarded Green's func- 
tion, which is determined by the operator □, (ii) the 
source term p, and (iii) the initial data on a Cauchy sur- 
face. Our goal is to represent with perturbation theory a 
physical situation which is as close as possible to that of 
the Misner data. For (ii), we take the source term p to be 
given by the time symmetric trajectory depicted in Fig. 
|l|, as discussed above. For (iii), we shall first study the 
case of no outgoing wave from the past horizon r = 2m 2 
and no incoming wave from past null infinity r = oo. We 
shall return to a more careful consideration of this point 
later. 

With (ii) and (iii) so specified, it is straightforward to 
write down the Zerilli function expressed for each multi- 
pole £ as the Fourier integral 

/>oo 

ip p (t-r*) =2y/2m 1 / du Z e (uj)e- l ^ t - r '\ (3.11) 



where 



^-1)1(1+1)^+2) 



2/ 



uA^uj) (£-l)t(£ + l)(£ + 2) 



12irri2UJ 



in 



ro 



2m 



2 J - 



drj 



yjl- 2m 2 /r + r 



1 — 2m2/r 



xe^HX^r). 



(3.12) 



Here, r\ is a parameter along the trajectory, such that 
r = (r /2)(l + cosry), r = dr/dr = -(2m 2 /r ) 1/2 (l + 
cos?7) -1 sinry, with r denoting proper time, and 



1 



-Jo 



2m 2 



4m 2 
ro 



r\ + sin rj 



-2rri2 In 



y / r /2m 2 — 1 + tan \r) 



\Jroj2m 2 — 1 — tan \r] 



(3.13) 



r\ increases monotonically along the trajectory, and takes 
the values TVh, where 7]h — cos~ 1 (4m 2 /r — 1), at 
the past and future horizons, respectively. The function 
Xi(oj;r) satisfies the Regge- Wheeler equation, 



d 2 
dr* 2 



2m 2 



1(1 - 1) 6m 2 



>X e = 0, 
(3.14) 

with boundary condition Xi(uj;r — > 2m 2 ) ~ e~~ lu}r . At 
large values of r, the Regge- Wheeler function becomes 
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X t (uj;r oo) ~ A t (w)e- iur " + 0(e^ r *), which de- 
fines the constant ^(w) appearing in Eq. ( 3.12 ). Finally, 
T((lu) is the differential operator 



r t (w) = 2(1 - 2m 2 /r + iwr)r(l - 2m 2 /r) 

+ (1 - 2m 2 /r) [£(£ + 1) - 6m 2 /r] 
+ 2iujr(l — 2m 2 /r + iwr). 



(3.15) 



In Fig. ^a we show the I = 2 component of -0 p for the 
case r = 15m 2 ; the Zcrilli function is plotted as a func- 
tion of coordinate time t, and is measured by a detector 
situated at radius r — 40m 2 . At t — — 29.5m 2 a sharp 
feature appears, corresponding to the particle emerging 
from the past horizon (t w in Fig. ||). The white hole 
subsequently goes into quasinormal ringing. In Fig. ^a 
the dotted line represents a pure quasinormal-mode sig- 
nal, the superposition of the first two (the least damped) 
1 = 2 quasinormal modes of a Schwarzschild spacetime 
of mass m 2 , with amplitude and phase determined by 
matching to the white-hole ringing. The dotted line was 
plotted starting from t = — 22m 2 , but it is practically 
indistinguishable from the solid line in the early part. 
The good agreement between the wavelengths confirms 
that the early-time portion of ip p is indeed quasinormal 
ringing. The ringing eventually stops, and the subse- 
quent low-frequency signal at 30m 2 < t < 80m 2 is the 
bremsstrahlung radiation coming directly from the par- 
ticle as it emerges from the strong field region. As the 
particle moves outside r = 3m 2 , which is the peak of 
the potential barrier for the scattering of gravitational 
waves, the bremsstrahlung radiation can reach the ob- 
server without much scattering by the spacetime curva- 
ture. The ray emitted at the turning point reaches the 
detector at time t — 27m 2 (cf. point t s in Fig. ||). 
The particle then turns around and falls back. As the 
particle once again goes through the strong field region 
while approaching the future horizon of the black hole, 
quasinormal modes are excited once more. It can be seen 
from the figure that the black-hole ringing has a smaller 
amplitude than the white hole's. 

In Fig. [!)), we show the t = 4 component of ip p 
for the same case. We see that the radiation is again 
consist of the same three components, except that the 
bremsstrahlung contribution to £ — 4 is much weaker. 
The quasinormal mode nature in the early part of the 
waveform is clear, with the dotted line giving the quasi- 
normal mode fit (first two £ = 4 modes) . The dotted line 
begins at t — — 25m 2 and is nearly indistinguishable from 
the actual waveform. Notice that the black hole ringing 
has a much lower amplitude. 

The waveforms given in Figs. ||ab point to the fact 
that the Misner data represent more than just two black 
holes in time symmetric trajectories. If there were no 
wave coming from the past null infinity in the Misner 
data, we expect that there would be "white hole ringing" 
as shown in Figs, ^ab, for black holes in time symmetric 
motion [El). As the Misner data produces no whitc- 



( = 2, r„ = 15m, 



white/black hole ringing 
white hole QNM fit 




-50.0 0.0 50.0 100.0 150.0 200.0 
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1 = 4, i„ = 15m, 



white/black hole ringing " 
white hole QNM fit 
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Time 
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FIG. 4. The I — 2 and I = 4 components of ip p for the case 
ro = 15rri2; the Zerilli function is plotted as a function of co- 
ordinate time t, and is measured by a detector situated at ra- 
dius r — 40m2 . There are three components in the waveform: 
white hole ringing, direct bremsstrahlung radiation, and black 
hole ringing. The dotted lines show the quasinormal mode fit 
to the white hole ringings. 
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Particle-Membrane Waveform for r„=15M 



hole ringing for t > 0, independent of the location of the 
detector, the initial data given at t = must contain the 
right amount of waves traveling inward to cancel the out- 
going "white hole ringing" . The time symmetry implies 
that there is no net flux at t = 0. The fact that there 
must be waves coming in from past null infinity is also 
guaranteed by time symmetry: there are as many waves 
coming in from past null infinity as going out to future 
null infinity. The same argument can also be applied to 
waves in and out of future and past horizons. 

These considerations suggest that, for a meaningful 
comparison to the Misner data, the perturbation calcu- 
lation we sh ould carry out is not exactly the ip p given 
by Eq. (3.11) above, which assumes no out going waves 
from the past horizon ( H~) and no in coming waves 
from past null infinity ( Instead, nontrivial initial 

data [component (iii) in the specification of the retarded 
integral, Eq. ( 3.1Cl| )] should be imposed on H~ and I~ to 
ensure that there is no net flux at t = 0. In particular, 
the Cauchy development of these data should cancel the 
white hole ringing. 

Obviously, the construction of such data on H~ and 
I~ is difficult, if at all practical. We circumvent this 
difficulty by directly subtracting the "offending" compo- 
nents of ipp shown in Figs. ||. [What this subtraction 
corresponds to at H~ and I~ can in principle be deter- 
mined by integrating backward in time. But that is not 
the concern of this paper.] The result of this subtraction 
is shown in figure |^. The components that we want to 
subtract can readily be identified: First, the white hole 
ringing. Second, we note that the portion of the wave 
given as a dashed line has retarded time less than 0. It is 
emitted by the particle as it is flying out before reaching 
ro (t < t s in Fig. |^), and therefore has no corresponding 
part in the Misner data. [As the scattering due to the 
potential is weak once the particle is outside the peak of 
the potential barrier at r — 3m2, this contribution to the 
waveform can be identified by its retarded time. This is 
not so for quasinormal ringing, which is a multiple scat- 
tering phenomenon.] 

These subtractions based on physical understanding of 
the system, although not expected to yield mathemati- 
cally exact time symmetric spacetimes, make it possi- 
ble to use perturbation theory to construct a waveform 
which is comparable to that computed numerically. Both 
its similarities and differences with the numerical results 
shed light on the physical meaning of the Misner data. 

The next step in the construction of the waveform 
is to extrapolate our results to the equal-mass case. 
One reasonable choice is to replace mi wi th \l (t he re- 
duced mass), while all m-2 in Eqs. (3.11)— (3.15) with 
M = mi +1TI2 — 2rri2 (the total mass). This respectively 
fixes the amplitude of the Zerilli function, and the units 
of the time coordinate t. We also correct for the internal 
dynamics of the "particle" (now imagined to be a black 
hole) by m ultiplying the Zerilli function by (Th) 1 ^ 2 [c.f. 
Eqs. ( |3.7[ )]. The result is our particle- membrane wave- 
form, -0pm, displayed in Fig. || for the I — 2 component. 




50.0 100.0 
Time 



FIG. 5. We show the i — 2 component of the parti- 
cle-membrane waveform for the case ro = Vbuii- The present 
numerical relativity code cannot evolve initial data sets with 
such large separations long enough to determine the full wave- 
form. 



This figure corresponds to the case where the maximum 
separation between the two black hole is ro = 15m2- 

It was useful to consider the case r = 15m 2 be- 
cause the gravitational-wave signal shows (before sub- 
traction) three relatively well separated portions: white- 
hole ringing, particle bremsstrahlung, and black-hole 
ringing. However, and unfortunately, numerical results 
are not available for such a large initial separation. (Dif- 
ficulties associated with the Cadez coordinate system 
make the numerical evolution unreliable when the initial 
separation is large). So Fig. || can be regarded as a "pre- 
diction" of what the numerical evolution of the Misner 
initial data should produce when pushed to such initial 
separations. 

To compare "0pm to existing numerical waveforms, we 
consider the cases ro = 6.8M (corresponding to a value 
fj,0 = 3.0 for the Misner parameter) and ro = 5AM 
(/io = 2.7). For such initial separations, the three por- 
tions of the waveforms are not cleanly separated. In 
Fig. ||a we consider the case fXQ = 3.0, and display the 
1 = 2 component of three distinct waveforms; the wave- 
forms are measured by a detector situated at r = 35M. 
The first i s ip v represented by the dotted line, obtained 
from Eq. (3.11) (with extrapolation mi = /i,m-2 = M 



). The second is tppm represented by the dashed line, 
obtained from ip p by subtracting the white-hole ringing 
and including the corrections due to internal dynamics. 
The third is the numerical waveform, V'num as solid line. 
In Fig. |p3 we do the same for the case fio = 2.7 (the 
detector is situated at r = 20M). 

In Fig. ^a, the white-hole ringing produces peaks at 
around t/M = (19,35,52); the last peak overlaps with 
the black-hole ringing. The bremsstrahlung radiation is 
also contained in the overlap. Because of the overlap 
with the bremsstrahlung radiation, the white-hole ring- 
ing cannot be subtracted efficiently, and the resulting 
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FIG. 6. Comparison of the I = 2 particle-membrane wave- 
forms with the numerical results for separations /io = 3.0 and 
fio — 2.7. The dotted line is the result of the particle tra- 
jectory before the "white hole" part has been subtracted, the 
dashed line is the particle-membrane waveform corrected for 
the white hole ringing, and the solid line is the full numerical 
result. 



-0.2 




-0.4 
-0.6 



-50.0 0.0 50.0 100.0 150.0 200.0 

Time 

I = 2, n = 2.0 




-0.2 - 



-0.4 i 
-0.6 - 

-50.0 0.0 50.0 100.0 150.0 200.0 

Time 

FIG. 7. Comparison of the 1 = 2 particle-membrane wave- 
forms with the numerical results for separations /io = 2.2 and 
2.0. 

"0pm remains contaminated. The contamination is less 
in Fig. ^b, because less bremsstrahlung is present in the 
signal (the magnitude of the bremsstrahlung radiation is 
determined by the speed of the particle in the part of 
the trajectory with r > 3M). Moreover, in this case the 
white-hole ringing is more in phase with the black hole's, 
hence the contamination less noticeable. We do not see 
any "hump" like the one that appeared at t = 52M in 
Fig. 6a. Instead, the contamination showed up as a 
slight lengthening of the wavelength of "0pm . Apart from 
these effects, the agreement between ip pm and V'num is 
reasonably good. 

The flat plateau at early times (t = to 37 for [i — 3.0 
and t = to 25 for // = 2.7) in lmm is a non-radiative 
component in the Misner data. It's magnitude depends 
on the location of the detector, chosen to be located at 
35 M for fi = 3.0 and 20 M for fx = 2.7, with the former 
one having a much smaller amplitude. Should the de- 
tector be put further out, this component would be even 
smaller. After the quasi-mode ringing sets in, we see that 
V^mm agrees very well with ippm both in their phases and 
amplitude. This is remarkable, in view of the fact that 
the matching involves no adjustable parameters. 
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In Figs. 0a and 0b we push the particle-membrane ap- 
proach to even smaller initial separations. In Fig. 0a we 
display the waveforms corresponding to the case /lo = 2.2 
(ro = 3.7 M), while in Fig. 0b we consider the case 
fio = 2.0 (ro = 3.2M). We sec that, for such initial sep- 
arations, the waveform ip p obtained from the perturba- 
tion calculations (dotted lines in Figs. 0a and 0b ) has its 
white-hole and black-hole ringing parts merged to a large 
extent. This is expected, because r is very close to 3M, 
the location of the peak of the potential barrier surround- 
ing the hole. For such cases, in order to separate out the 
white hole ringing, we match ip p to two sets of quasinor- 
mal modes which arc excited at two different times. The 
first set is denoted V\v ("w" stands for "white"), and is 
composed of the first two £ = 2 quasinormal modes. The 
second set which sets in at a later time is denoted iph 
(with "b" standing for "black"), has the same frequen- 
cies. With the amplitude of ip w determined through such 
a matching, we obtain ip pm by subtracting ifj w from tp p . 
In Figs. 0a and 0b, V'pm is represented as dashed lines. 
The agreement between ip pm and ~ ^Anum is quite satis- 
factory except for the initial pulse. The disagreement in 
the initial pulse is due to the difficulty in determining ip w 
from the merged white hole-black hole ringing. For the 
case of no — 2.0, we find that the the white hole ringing 
■0w, which is to be subtracted away, is quite sensitive to 
the fine details in the form of ip p . This signals the be- 
ginning of the breakdown of the method. As pointed out 
above, we do not expect our particle-membrane approach 
to be applicable for small initial separation between the 
two holes. For initial data with /io = 1.8 or smaller, the 
two holes are surrounded by a common event horizon. 
This breaking down of the method at around /io = 2.0 is 
consistent with the results on energy radiated as shown 
in Fig. The particle-membrane result also begins to 
deviate from the numerical result at around /io = 2.0 
there. 

In Figs. 0, we compare the £ = 4 components of the 
waveforms for various cases. The dotted lines give ip p . 
The dashed lines give the particle-membrane waveform 
"0pm extracted from ip p . They are superimposed on the 
numerical results (solid lines). Fig. 0a is for the case 
^o = 3.0, while Fig. 0b is for pt = 2.7. We see that 
the white-hole ringing has much larger amplitude than 
the black-hole ringing (which are shown blown up in the 
insets), making the extraction more difficult. However, 
we believe that the pm so obtained is still more reliable 
than the numerical ones, which are considerably larger 
in amplitude. In the next section we shall see a similar 
situation for small /i values under the "close approxima- 
tion" . Although the numerical calculations give wave- 
forms agreeing very accurately with the semi-analytic 
methods for the I = 2 components, the agreement for 
the £ = 4 components is much less satisfactory. 

We conclude this subsection with the following re- 
marks. In view of the complications associated with the 
removal of the white-hole ringing from the perturbative 
waveform, it is tempting to rule that only the infall part 
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FIG. 8. Comparison of the £ = 4 particle-membrane wave- 
forms with the numerical results for separations = 2.7, 3.0. 
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of the trajectory should be included as source from the 
first place. Thus, one might ignore the first half of the 
motion, in which the particle travels from the past hori- 
zon to the t urnin g point at r = tq. One might then 
modify Eq. ( [3.12 ) so that the integral is now evaluated 
between 77 = and r) = T}h- However, such a truncation 
would produce nonsense: It corresponds to the sudden 
creation of a particle at r = 7*0. Such a violation of 
the Einstein's equations, which guarantee the conserva- 
tion of energy- momentum, produces unacceptable results 
even at the linearized level considered here. This has also 
been observed previously in Ref. |24| . 

We sh ould n ote th at our present calculation, as given 
by Eqs. ( 3.11 )^( 3.15 ), has also violated the Einstein equa- 
tions. The Teukolsky equation was integrated assuming 
no outgoing radiation from the entire past horizon. How- 
ever, this is not possible with a particle flying out of it 
(c/. ; Fig. 3). If one imposes a no out going wave condi- 
tion as initial data for the portion of the horizon at times 
later than the emergence point of the particle (point E in 
Fig. 3), then we cannot impose the same on the portion of 
the horizon before point E. Otherwise, Einstein's equa- 
tions are violated at the point E, with the source term 
p being a 5 function having support there. To satisfy 
Einstein's equations, one has two options. One can let 
the particle trajectory extend all the way back to the ini- 
tial singularity, eliminating the need of specifying bound- 
ary conditions on the past horizon. While allowing the 
passage through the coordinate singularity at r = 2m,2 
should be straightforward, the difficulty is that there is 
no a priori clear way to pick the initial data at the past 
singularity that best corresponds to the Misner data. An- 
other option is, while imposing a no out going radiation 
condition for the portion of the horizon at times later 
than the emergence point E of the particle, we integrate 
the Teukolsky equation backward along the past horizon 
across the point E to determine the suitable data on the 
earlier part of the horizon. However, we note that do- 
ing this extra work to get a more "correct" set of initial 
data will not affect much our final waveform ip pm . This 
is because for the data so determined on the early part 
of the past horizon, the high frequency components of it 
will propagate out unscattered, reaching the detector at 
a very early retarded time and hence cannot not interfere 
with the waveform we want to extract, while the low fre- 
quency components of it will be trapped by the potential 
barrier and appear later as quasinormal ringing of the 
white hole, which we subtract away anyway. There will 
only be a negligibly small part multi-scattered by the po- 
tential outside and away from its peak (r > 3 M) that 
can reach the detector at the time of our interest. Hence 
we see that, for our present purpose, carrying out this ex- 
tra step of imposing correct data on the early part of the 
past horizon can have only negligible effect on our final 
waveform tp pm . This also explains why violating the Ein- 
stein equations in our present setup is not troublesome 
(whereas it is troublesome if the trajectory is truncated 
at the turning point), as evident in the close matching of 



our results to the numerical results shown above. 

To sum up, in this section we have calculated the wave- 
form generated by a particle in time symmetric motion, 
with no wave coming from the past horizon and the past 
null infinity. We then subtract away the part of the ra- 
diation emitted by the particle before the time symmet- 
ric point t — 0. The resulting waveform, with a cor- 
rection factor (of order unity) put in to correct for the 
internal dynamics of a black hole, is compared to the 
waveform obtained by numerically evolving the Misner 
data. We find that for [i > 2.0, the agreements in the 
phases, frequencies and amplitudes of the two waveforms 
are satisfactory. This implies that the Misner data repre- 
sents, to a good approximation, such a physical situation, 
namely, two throats in time symmetric motion with ra- 
diation from the t < part of the trajectories balanced 
by waves from the past horizon and the past null infinity. 
For cases with smaller mu values, we have difficulty in 
identifying, and hence subtracting the radiation emitted 
by the "particle" before the time symmetric point, and 
we cannot obtain a waveform from this semi-analytic ap- 
proach. Fortunately, for smaller mu values, we have de- 
veloped a different semi-analytic treatment which we now 
turn to. 



IV. PERTURBATION THEORY FOR THE CLOSE 
LIMIT 

A. Formulation 

When black holes start sufficiently close to each other 
the analysis of the radiation generation can be consider- 
ably simplified. From numerical evolution computations 
it is known that for fiQ less than ~ 1.8 a common horizon 
initially surrounds both throats of the Misner geometry 
]l8t . We describe here how the spacetime exterior to this 
horizon can be viewed as a distorted Schwarzschild geom- 
etry and can be treated with methods from perturbation 
theory. 

Our starting point for the perturbation analysis is the 
Misner geometry |J, 

^Misner = aVtlisnor + <V + SU1 2 T] d(j) 2 ] , (4.1) 

where 



fMii 



E 



(4.2) 



n— — oo 



\J cosh(/i + 2nfio) — cos r\ 
and where the coordinates have the range 

-MO <M<M0 < 77 < 7T < </> < 2?T . (4.3) 

This geometry represents an asymptotically flat three ge- 
ometry with ADM mass 



M = 4a V — = 4a£i 

z — ' sinhnun 



(4.4) 
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The geometry has two "throats" which end at fi = ±/io- 
These throats can be considered to be joined to form a 
wormhole, as in the original Misner paper, in which case 
the Misner solution is viewed as periodic in fi, with period 
2/io- Alternatively both throats can be extended into a 
second asymptotically flat space [p5| , in which case the 
period is 4/_to and range /xo < /x < 3/xo may be considered 
to describe the mirror image solution with the second 
asymptotically flat space. 

The nature of the Misner metric as a single perturbed 
hole becomes clearer if new coordinates R, 9 are intro- 
duced as if [i, 77 were bispherical coordinates being trans- 
formed to spherical polars: 



R = a\ 



I cosh fi + cos 77 
cosh /1 — cos 77 



tan 6* 



smr7 
sinh [i 



(4.5) 



If we want to restrict our R, 9 coordinates to the region 
— Mo _ M _ Mo we must in principle remove the interiors 
of two "circles" corresponding to |/x| > [j,q. 

In terms of these new coordinates, the Misner geometry 
takes the form 



ds 



Mi: 



$ 4 (R, 9; fj, ) (dR 2 + R 2 [d9 2 + sin 2 9 d(j> 2 } ) 



(4.6) 



with the conformal factor $ given by, 



(4.7) 



-1/2 



(4.8) 



$ = 1 + 2_, V cosn m ~ cos v x 

n/0 

{2 sinh 2 n/io + sinh /x sinh 2n/io 
+(cosh fj, — 1) cosh 2nfiQ + 1 — cos rj} 

This can be rewritten as, 

$ = 1 + 5 ^ [(1 + <5 2 ) sinh 2 n/j 

n/0 

+5 cos 9 sinh 2n/j + 5 2 ] 

where S = a/R = M/4KS±, and where we have elimi- 
nated a by using (PI). 

The sum in ( |4. 8 ) is recognized as the generating func- 
tion for the Legendre polynomials, so the conformal fac- 
tor can be rewritten as 

$ = 1 + 2 Ke( f io)(M/R) e+1 Pi (cos 0) . (4.9) 

£=0,2,4... 

The only /x dependence occurs in the Kg coefficients, 

. . 1 v—* (cothn/xo) 

K ^o) = — t+1 }^ 



-J sinh n/xg 



(4.10) 



and this dependence is shown in figure |9| 

If the I = term in the sum for $ is explicitly evaluated 
(note that ko(mo) — 1/4) the result is 



1=2 
1=4 



FIG. 9. The K( coefficient as a function of /xo- 



with 
T= 1 + 2(1 



M 
2R 



M 
2R 



(4.11) 



KM(M/R) e+1 P e (cos9). 



1=1 A.. 



(4.12) 



Since R is analogous to the radial coordinate in 
isotropic coordinates for a spherical geometry, we in- 
troduce a Schwarzschild-like radial coordinate r by the 
transformation 



R 



(4.13) 



that relates the two coordinates in the Schwarzschild ge- 
ometry. With this transformation we arrive, finally, at 



^ s Misncr — ? ( r > 



dr 2 



1 - 2M/r 



l dn- 



(4.14) 



in w hich t he r and dependence of T are given by ( 4.12| ) 
and (4.13). It is easy to show that a choice of time coor- 
dinate t can be made so that the 4-geometry generated 
by the Misner initial data takes the form 



ds 2 



2M 
r 



+T{r,9f 



dt 2 
dr 2 



(4.15) 



1 - 2M/r 



r 2 dn 2 



at t = with all first time derivatives of the metric van- 
ishing. 
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The function T may therefore be thought of as con- 
taining the mathematical description of how the Misner 
geometry initially deviates from the Schwarzschild geom- 
etry. 



In treating (4.15) as a distorted Schwarzschild geome- 
try we are takin g the range of r to be 2M < r < oo and 
hence, by ( 4.13| ), our range of coordinates must include 
all points with R > M/2. But according to the transfor- 
mation in ( [Ow the values of R at r\ = are limited (by 
the condition \fi\ < fio) to values for which 



R > acoth/xo/2 = (M/4£i) coth/i /2 



(4.16) 



Thus our coordinates can reach down to R — M/2, while 
corresponding to > fio only if coth/io/2 < 2Si, and 
this condition turns out to be satisfied only for [j,q < 
1.191. We ignore this restriction and find (see below) 
that the perturbation approach gives reasonably accu- 
rate answers for values of hq somewhat larger than 1.191. 
The explanation is that all the strong non-spherical de- 
viations are very close to the horizon and end up not 
having important consequences for emitted radiation. 

Note that the Misner geometry with the Kg set to zero 
would be precisely at = const, slice of the Schwarzschild 
geometry. The deviatio ns fro m Schwarzschild are de- 
termined by the k^'s in ( 4.12 ). If they are sufficiently 



small then the Misner geometry may be considered ini- 
tially to be a perturbed Schwarzschild geometry. Since 
the Schwarzschild geometry is stable, sufficiently small 
initial perturbations will remain perturbations, and the 
spacetime generated from the Misner initial data will be 
a perturbed Schwarzschild spacetime. A crucial point is 
that the form of the initial data interior to the initial 
horizon cannot affect the evolution of the exterior. We 
can, therefore, view the generation of outgoing waves — 
a process confined to the exterior of the horizon — as 
nearly spherically symmetric if the initial data is nearly 
spherically symmetric only outside the horizon. More 
specifically: if M/R is of order unity, or smaller, outside 
the initial horizon, and if the Ke(no) coefficients are small, 
then the evolved spacetime outside the horizon will be a 
perturbed Schwarzschild geometry. 

This picture makes sense only if the coefficients Kt's 
become small when fi$ — > 0, and if this limit corresponds 
to the "close limit" in which the separation of the two 
black holes vanishes. To show the former we notice that 
(cothn/zo)V sinhn/^o « (nfj,o) e+1 for n/i << 1. We as- 
sume that this approximation is valid for small fj,o, and 
we keep terms in summations only up to N ~ 1/mo- With 
these approximations we have, for I > 1. 



(cothn/io) 
^— ' sinh nun 

n—l 

and similarly 



(iM 



1 N 1 



C(*+i) 



Mo 1 



(4.17) 



sinh n^o Mo ~ n 



1 N 1 



In TV 
Mo 



I In Mo | 

Mo 



These approximations in (4.10) give us 



k^(Mo) 



|41n Mo r 



(4.18) 



(4.19) 



Although the steps used to derive this result were very 
rough approximations, the result is in good agreement 
with numerical values of the Kg(fio) s. 

We turn now to the question of the separation of the 
throats in the limit of small /io- The proper distance L 
between the throats (more specifically, the distance from 
fx = — /iQ to +/io along the line rj = tt) can be written as 
the sum [Ef| 



L = 2a 1 + 2/x 



sinhn/j,o 



(4.20) 



Here it will be more useful to use an equivalent result in 
which L(po) is given implicitly in terms of K and E, the 
complete elliptic integrals of the first and second kind: 



Mo 



{4,a/Tr)K \{l - fc 2 ) 1/2 l E{k) 



irK(k)/K Ul-k 2 ) 1 ' 2 



(4.21) 



(4.22) 



If /io is small, then this forces k <C 1 so that E(k) 
K(k) re tt/2, and 



L=(7r 2 a/Mo) [l + 0{k 2 )} 
with k ps 4exp (— it 2 /2{i ). We then have 



L 

M 



4/i Si 4|ln/i | 



(4.23) 



(4.24) 



and hence the separation decreases as Mo decreases. 

We can then consider e = 1/| In /xq| to be a pertur- 
bation parameter, with which we can evolve the Misner 
initial using perturbation theory. The form of the se- 
ries in ( 4.12| ) appears to give a multipole decomposition 
of the deviations from spherical symmetry. This is not 
quite true. The metric perturbations involve T raised 
to the fourth power. When this is done terms of differ- 
ent multipolarity in T mix to give each multipole of T . 
The quadrupole term in J 74 , for example, will contain 
contributions from the product of the £ = 2 and £ = 4 
terms in T , from the square of the £ = 2 term, etc. The 
total quadrupole term will contain products of the form 
K ii K8K2K4, and so forth. This complexity causes 
no difficulty in practice. Of all the contributions to the 
quadrupole in JF 4 , that of lowest order is the 0(e 3 ) term 
linear in k 2 ; the next lowest order contribution C(e 6 ) is 
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that from k\ . If we are evolving the initial data with the 
evolution equations of linear perturbation theory for the 
quadrupole then we must keep only the part of T A that is 
linear in k-x. It would be inconsistent to keep any higher 
order terms since we are ignoring the contributions of 
order C(e 6 ) due to the nonlinear evolution of the initial 
data. A variation of the same argument applies in the 
case of other multipoles. If we are using the £-pole linear 
equations to evolve the £-pole initial data, we must keep 
only the term in JF 4 that is linear in ki- All other £-pole 
contributions will be higher order in e. Thus if we are 
only interested, for each ^-pole order, in the terms that 
can be treated with linearized perturbation theory, we 
may write T as 



1 + 8(1 + Kt(j*o)(M/R) i+1 P t (coB9), 

(4.25) 



£=2,4.. 



where 



Pe(cos ( 



4tt 



2^+1 



-Yi 



an 



(4.26) 



To evolve these even parity perturbations we use the 
formalism and the notation of Moncrief |l6f| . From 
Eqs. (5.1)-(5.8) of that pape r, appl ied to the perturbed 
spacetime of our Eqs. ( 4.15| )-( [4.25 ), we find that Mon- 
crief's perturbation functions h\ and G vanish, and 



H 9 



K = g e (r; [i Q ) = 8 



Air 
2£+l 



(M/R) 



tt+i 



M/2R 
(4.27) 



Qt-\i 



and Moncrief 's q\ is 



2M 

qi = 2r 1 



+£(£ + \)rQ t . 
We then define, for any £ 

(e + 2)(£-i) 




pert 



(£ + 2)(£-l) + 6M/r 



(4.28) 



(4.29) 



For 1 = 2 note that V'pert is identical to the ip of [£7j . Note 
also that the definitions here for th e m et ric p erturbations 
are closely related to those in Eqs. (2^)- (2^4), except that 
the normalization used for ip is different (see below ). 

A t thi s point it should be observed that in Eqs. (4.27) 
to (p|) the metric perturbations are proportional to 



Ki(^o) and that there is no other /io dependence in the 
perturbations. We may therefore view these expressions, 
for each £-pole moment, as first-order perturbation the- 
ory in K^(/i ). From a formal point of view we may con- 
sider that we have a family of spacetimes parameterized 
by Ho, and we are approximating the metric by 



g a f3 ~ ffa/3| K( ,( po ) = o 



+ 



«f («i)=0 



(4.30) 



It is important to realize that we could just as well do 
perturbation theory using another expansion parameter 
Ke(p-o) that agrees to linear order with ^^(/io); we might 
for example use Ki(^o) = nei^o) + [nei^o] ■ I n terms of 
this parameter our first-order approximation would be 



g a /3 ~ 9af3\ ke ( Mo ) =a 



+ 



X Ki(po) 



«<?(a<o)=0 



dg a /: 



x k e (p ) (4.31) 



The perturbations here differ from those in (4.30) by 
the factor ki (no)/ni(p-o). The approximation given by 
first-order perturbation theory, then, is dependent on our 
choice of expansion parameter. This, of course, is a re- 
flection of the fact that first-order perturbation theory 
must be uncertain to second-order in the expansion fac- 
tor, but it should be kept in mind that the remarkably 
good predictions of perturbation theory at unexpectedly 
large values of /iq are to some extent due to the particu- 
lar (though rather natural) use of Ki(po) as the effective 
expansion parameter. 

The function ip per: t is evolved with the Zerilli equation 



pert 



d 2 ip 



pert 



dt 2 



Or 



*2 



K(r*)V : 



pert 



= 



(4.32) 



in which r* is a "tortoise" coordinate such that r* — > oo 
corresponds to spatial infinity and r* —oo corresponds 
to the horizon, 



r* = r + 2M\n (— 1 

\2M 

and the Zerilli potential is given by, 



(4.33) 
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(4.34) 



and A = l + 3M/2r. 

The initial form of V'portj given by Eqs. ( 4.27 )-( 4.2E ), is 



shown in Fig. |10| along with the potential for the Zerilli 
equation. Both the and £ = 2 and £ = 4 figures are 
shown. 

It is straightforward to show that, in terms of ippert, 
the radiation power is given by 



dPower 
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/ <#pcrt \ 2 
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and therefore, 



Power = 
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FIG. 10. We show the initial perturbation waveform, com- 
puted at t = directly from the Misner initial data, and the 
Zerilli potential for I — 2,4. 
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167T 



#pei 

dt 



(1-2)! 
{t + 2)\ 



(4.36) 



Since the waveforms and power depend on /io only 
through factors of Kg(/J,o), we can compute with set 
to unity, and get the correct results for any /io by multi- 
plying waveforms and power expressions respectively by 
Ke(no) and n e (fi ) 2 . 

The above formalism has been used to generate the 
waveforms and the energies that are presented below and 
compared to the results of numerical relativity. In the 
comparisons we must take account of the fact that dif- 
ferent normalization conventions than those above have 
been used to define the wave function ip num . The compar- 
isons will be made by converting the perturbation wave- 
form according to 



Ij'n 




pert 



(4.37) 



B. Comparison of results 

Probably the most important comparison to be made 
is the radiated energy computed by the different methods 
for dealing with the outgoing radiation. In Fig. [ll] this 
comparison is given. The numerical relativity supercom- 
puter results are shown along with error bars indicating 
the range of energies found by extracting waveforms at 
different radii. There are no analogous formal errors for 
the close-limit or the perturbation-paradigm approach. 

The results show remarkable agreement (already noted 
in p7| ) between the close- limit prediction and the re- 
sults of fully nonlinear numerical relativity. This agree- 
ment is reasonably good even up to /io = 2, where the 
energy results differ by only 15%. Above /xq = 2 the re- 
sults quickly diverge, with the close-limit results seriously 
overestimating the radiated energy. For these cases, how- 
ever, the perturbation-paradigm method gives excellent 
agreement with computed results (as has already been 
noted in j5|||). 

The energies shown in Fig. [ll] are strongly dominated 
by the quadrupole contribution. It is instructive to look 
at the £ = 4 energy predictions of numerical relativity 
and of the close limit. In the case of the close limit, the 
£ = 4 energy is computed in the manner explained in 
Sec. IVA. In the numerical relativity results the £ = 4 
results must be extracted by fitting the angular pattern 
as explained in Sec. II. The numerical £ — 4 results lack 
the clear trend, seen in Fig. [jj], of agreement at small /io, 
monotonically growing worse with increasing /zo- 

A better understanding of the sources of disagreement 
of the method follows from an examination of the "wave- 
forms," the time dependence of the perturbation func- 
tion ip at a constant radius. Figure O shows a series 
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— close limit 
- pert-para 
■ numerical rel. 



FIG. 11. Comparison of the close approximation energy 
with the full numerical results and the particle-membrane ap- 
proximation 



L=4 close limit 
L=4 numerical rel. 
L=2 close limit 
L=2 numerical rel. 



FIG. 12. Comparison of energies for both the t = 2 and 
I = 4 modes. The numerical relativity results are denoted by 
symbols showing the specific case computed, while the close 
approximation results are denoted by lines. 



of waveforms "observed" at r = 25M (where here and 
below M refers to the ADM mass of the initial data), 
and show the astonishingly good agreement at small (j,q, 
agreeing in form not only in the region dominated by 
"ringing" at the quasinormal frequency, but also agree- 
ing at early times. As fj,o increases, the curves continue to 
agree in general character, but disagree in the amplitude 
of the quasinormal ringing. Only a single interesting fea- 
ture (aside from the agreement!) appears in these curves. 
There is a consistent phase drift between the close-limit 
waveforms, and the numerical relativity waveforms; the 
late-time quasinormal ringing of the numerical relativity 
waveforms has a frequency too low by 10-20%. (The cor- 
rect values of the quasinormal frequencies are well known 
from other calculations. See e.g., |28f] ). As pointed out 
in Sec. II, this suggests that there is a systematic effect 
in the present numerical relativity computations causing 
this phase drift and leading to underestimates of radiated 
energy. Any increase in the numerical computed energy 
at large /xo will, of course, improve the agreement with 
the close-limit estimates. 

The comparison of £ — 4 waveforms, at r = 25M, 
shown in Fig. O, clarifies the trends in Fig. 12. For small 
/io numerical errors in the multipole extraction scheme 
result in waveforms which are clearly in error. Not only 
do the waveforms lack the expected quasinormal ringing, 
they have a non-physical trend at late times. These er- 
rors are so large simply because at small fj,o the radiation 
is overwhelmingly dominated by the quadrupole part. At 
/Jo = 1-0, for example, the i = 4 contribution is only 1% 
of the amplitude of the wave. As pointed out in Sec. II, 
extraction of this very small part is very sensitive to nu- 
merical noise. As /j,q increases the relative size of the 
i = 4 contribution increases, and the numerical error in- 
volved in extracting it decreases. For fio w 1.2 — 1.3 the 
numerical errors are small enough so that the waveforms 
show overall reasonable agreement (along with curious 
features at early and late times). At larger values of /iq 
it would be expected that the extracted waveforms would 
continue to be more accurate, but the lack of a monotonic 
increase in energy shows that the errors in the waveforms 
are still very large. Though the large errors make con- 
clusions uncertain, the results suggest that the linearized 
approach in the close-limit method has a smaller range of 
validity for I — 4 than for I = 2. Despite this somewhat 
smaller range of validity the point should not be missed 
that at the present state of the art in numerical relativ- 
ity, the close-limit waveforms and energies, for a range of 
fio, are the only reliable estimates available for I > 2. 

The waveforms presented above have been "observed" 
at radius r — 25 M. That is, the figures showed ip(r = 
25 M, t). To explore the sensitivity of the comparisons to 
different observation radii, in Figs, [lj we present fio = 1.2 
waveforms observed at different radii. Close-limit wave- 
forms are shown in Fig. p"5|a, and those of numerical rel- 
ativity in Fig. [l5|b. 

In these figures the phases of all waveforms have been 
made compatible with each other, and with Figs. O and 
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FIG. 13. Comparison of the £ — 2 waveforms obtained by 
the close approximation and those obtained by full numerical 
relativity for several cases. 



FIG. 14. Comparison of the waveforms for I = 4 from the 
close approximation and from numerical relativity. 
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r oSS =15M 

r„ h ,=35M 




0.0 20.0 40.0 60.0 80.0 100.0 

t/M 



r obs =15M 
r nh ,=35M 




FIG. 15. (a,b) Comparison of the waveforms "observed" at 
different radii for the close approximation and full numerical 
cases for /io = 1-2. (a) shows the results of the close approx- 
imation and (b) shows the full numerical relativity results. 
The phases have been adjusted for comparison of the wave- 
form shape at the different radii. 




FIG. 16. Comparison of the waveforms for the parti- 
cle-membrane approximation, numerical relativity and the 
close approximation for the case fj,o = 2.0. These three calcu- 
lations, containing no adjustable parameters, agree remark- 
ably well, although for this separation the close approximation 
is near the limit of its applicability. 

p"4| . Each waveform has been shifted by a time equal to 
the difference in observation values of r*, and the value 
of r* equivalent to r = 25 M. In Figs. |l5| we see that 
these phase corrected waveforms agree very well at late 
times, but have somewhat different initial shapes. At 
early times the waveform observed at smaller radius is 
larger than if observed at larger radius. This, of course, 
is a manifestation of the fact that at very early times 
what we are seeing is essentially the initial data, which 
unlike the outgoing radiation - falls off in radius. This ef- 
fect, the presence of a non-radiative part of the waveform, 
is significant up to around the first peak of quasinormal 
ringing. It is important to note that the differences in 
waveforms observed at different radii are present in both 
the close- limit and the numerical relativity results. In 
fact, the differences in the waveforms for different radii 
are much larger than the differences between the wave- 
forms computed by numerical relativity and by the close- 
limit approximation. 

The waveform comparison to this point has been be- 
tween the computations of numerical relativity and the 
close-limit approximation. We now tie this work together 
with the particle-membrane calculation described in sec- 
Although the two approximations are based 
as we have seen they nearly overlap 

In Fig. M 



on different limits, 

if the two holes are not too far or too close, 
we show results from all three approaches we have used 
in this paper for the case no — 2.0. The solid line shows 
the close approximation, the dashed line shows the re- 
sults from the full numerical relativity calculation, and 
the dotted line shows the particle-membrane calculation 
for the time symmetric particle perturbation. It is impor- 
tant to emphasize that the three graphs shown contain 
no adjustable parameters. 

The waveforms show good agreement at late retarded 
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time, when the behavior is dominated by quasinormal 
ringing. The waveform computed from the close limit 
has a higher amplitude, as this approximation is not as 
accurate for holes at such a large separation parameter. 
There is, however, a more significant difference in the 
shape of the early waveforms. Here it is the particle- 
membrane approach that d iffers from the other calcula- 
tions. As discussed in Sec. IIIB , for this fio — 2.0 case, 
the extraction of a waveform is difficult and is sensitive to 
the details of the extraction procedure. Overall, in view 
of the very different nature of the three calculations, and 
their different regimes of validity, the agreement among 
these results is remarkable. 



V. DISCUSSION AND CONCLUSION 

Conventional wisdom holds that events involving col- 
liding black holes are so nonlinear that they can only be 
studied using numerical relativity. We have shown in this 
paper that this is not necessarily the case. In particular, 
we have shown that the head-on collision of two black 
holes (whose starting point is the Misner initial data) 
is amenable to semi-analytic treatments giving reliable 
results in a considerable range of the parameter space. 

In the Misner data, the initial separation of the two 
black holes is determined by the value of the parameter 
/Jo- At present, our numerical techniques are capable of 
evolving the Misner data, and of extracting waveforms, 
for the range 0.7 < /io < 3.0, corresponding to initial 
proper separations L ranging from 1.51M to 7.92M. For 
Ho < 1.8 (L < 3.38M, corresponding to r = 2.82M), the 
black holes are surrounded by a common event horizon. 
In such cases the perturbative treatment ("close limit") 
described in Sec. IV provides a reliable framework for ap- 
proximating the true evolution of the spacetime, and the 
resulting waveforms match nicely those obtained numeri- 
cally. For fio > 1.8, the black holes are truly distinct, and 
we find that the waveforms can be well reproduced using 
the "particle-membrane" treatment described in Sec. III. 
This approach can also be used to "predict" the wave- 
forms produced during large- /io collisions, cases which 
should be amenable to numerical calculation in the near 
future. 

We note that in Sec. Ill B and IV B, the comparisons 
between the numerical and semi-analytic waveforms are 
made without using a single adjustable parameter. For 
the I = 2 component of the waveforms (the dominant 
component), the agreement between the numerical and 
semi-analytic approaches is remarkable. For £ — 4, ex- 
traction of the waveforms from the numerical data is 
less accurate, especially for smaller values of /io- As a 
result, and as shown in Sec. IV B, the numerical and 
semi-analytic £ = 4 waveforms differ significantly when 
Mo < 1-0- For larger values of /^o, the £ = 4 waveforms 
agree well in wavelength and in phase, but less so in am- 
plitude. This is illustrated in Figs. fL4| To explain this, 



we note that the waveform amplitude is more sensitive 
than other waveform attributes (such as wavelength) to 
the choice of parameters (such as resolution) used in the 
numerical calculation. 

What have we learned about the Misner data through 
the semi-analytic studies? For gravitational waves ob- 
served at large radii, the small /io cases of the Misner 
data actually represent just one single black hole with 
non-spherical perturbations, although the spacetime in 
the near field region can be very different. For larger 
values of /io, we have seen that the initial data sets rep- 
resent not just two throats in time symmetric motions. 
The waves coming from the past null infinity and the 
past horizon play important roles, namely, they cancel, 
to a large extent, the radiation emitted before the time 
symmetric point of the trajectories. We have also seen 
that the internal dynamics of the black holes do not have 
much effect on energies or waveforms observed at large 
radii. 

We are currently extending the semi-analytic ap- 
proaches to other types of black-hole events. We feel 
that such alternative approaches are an important com- 
plement to the direct numerical integration of the Ein- 
stein equations. Although these techniques cannot re- 
place numerical relativity, which will be the only means 
of computing detailed waveforms in more complicated 
spacetimes such as the full 3D inspiral of two black holes, 
they can augment it in a number of ways. First, they can 
be used to test the accuracy of the numerical results, as 
was illustrated in this paper. Second, they can predict 
results well in advance of the full-blown numerical treat- 
ment. Third, they are inexpensive computationally, and 
can therefore be used to search the parameter space (ini- 
tial separation, angular momenta, mass ratio, etc.) for 
potentially interesting phenomena, to be further inves- 
tigated using numerical relativity. Fourth, and perhaps 
most important, alternative approaches help provide a 
physical understanding of the numerical results. 
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